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Abstract: We investigate the short distance behavior of nucleon-nucleon (NN) potentials 
defined through Bethe-Salpeter wave functions, by perturbatively calculating anomalous 
dimensions of 6-quark operators in QCD. Thanks to the asymptotic freedom of QCD, 1- 
loop computations give certain exact results for the potentials in the zero distance limit. 
In particular the functional form of the S-state central NN potential at short distance 
r is predicted to be a little weaker than r~ 2 . On the other hand, due to the intriguing 
character of the anomalous dimension spectrum, perturbative considerations alone can not 
determine whether this potential is repulsive or attractive at short distances. A crude 
estimation suggests that the force at short distance is repulsive, as found numerically in 
lattice QCD. A similar behavior is found for the tensor potential. 
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1. Introduction 

In a recent paper Q a proposal has been made to study nucleon-nucleon (NN) potentials 
from a first principle QCD approach. In this field theoretic framework, potentials are 
obtained through the Schrodinger operator applied to Bethe-Salpeter (BS) wave functions 
defined by 

<Pe(t) = (0\N(x + r,t)N(x,t)\2N,E), (1.1) 

where |2N, E) is a QCD eigenstate with energy E (suppressing here other quantum num- 
bers), and N is a nucleon interpolating operator made of 3 quarks such as N(x) = 
e abc q a (^ x )q b (x)q c {x). Such wave functions have been measured through numerical simu- 
lations of the lattice regularized theory [[[], |2|, |3|, |4|. Although many conceptual ques- 
tions remain to be resolved, the corresponding potentials indeed qualitatively resemble 
phenomenological NN potentials which are widely used in nuclear physics. The force at 
medium to long distance (r > 2 fm) is shown to be attractive. This feature has long well 
been understood in terms of pion and other heavier meson exchanges. At short distance, 
a characteristic repulsive core is reproduced by the lattice QCD simulation jl]]. No simple 
theoretical explanation, however, exists so far for the origin of the repulsive core. For an 
approach based on string theories, see ref. ||. 
By writing 

(0\N(x + r,t)N(x,t) =J2 f ^(2N,mr,E\f n {r,E), (1.2) 

n=0 •* 

where |2N, mr,E) is a state with the energy E containing two nucleons and n pions 
(and/or nucleon-antinucleon pairs), we see that <£>£:(r) = fo(f, E). (Our normalization 
is (2N, nir, E\2N, n'-rr, E') = 2E5 nn /5(E — E').) We may thus interpret the wave function 
</?e(r) as an amplitude to find the QCD eigenstate |2N, E) in N(x + r, t)N(x, t)\0). 

The behavior of the wave functions <PE(r) a t short distances (r = \r\) are encoded in the 
operator product expansion (OPE) of N(x + r, t)N(x, t). An OPE analysis || of BS wave 
functions in the case of a toy model, the Ising field theory in 2-dimensions, successfully 
described the analytically known behavior. In this case the limiting short distance behavior 
of the potential does not depend on the energy (rapidity) of the state, and further it only 
mildly depends on energy (for low energies) at distances of the order of the Compton wave 
length of the particles. 

In this report we perform an operator product expansion (OPE) analysis of NN BS 
wave functions in QCD, with the aim to theoretically better understand the repulsive core 
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of the NN potential, (at least that of the measured BS potential). Thanks to the property 
of asymptotic freedom of QCD the form of leading short distance behavior of the coemcent 
functions can be computed using perturbation theory. A short summary of our results has 
been published in ref. Q . 

In sect. 2 we start with some general considerations on BS potentials, and sect. 3 
presents some standard renormalization group considerations. The anomalous dimensions 
of 3- and 6-quark operators are computed in sect. 4. Finally in sects. 5, 6 we discuss the 
application of the results to NN potentials. In appendix C we make a similar analysis for 
the 1 = 2 two pion system. For the convenience of the reader we give a brief summary of 
our results here. The OPE analysis shows that the NN central potential at short distance 
behaves as 

v„» - c j-^f'- 1 (1.3) 

for the S-state (L = 0) , where S and / are total spin and isospin of the NN system, 
respectively, and (3 SI is negative and explicitly obtained as 

(3 01 = , /3 10 = , (1.4) 

33 - 2JV f 33 - 2JVf 

where iVf is the number of quark flavors, and the overall coefficient Ce depends on the 
energy E. Unfortunately the OPE analysis is not as conclusive as that in the toy model 
referred to above, in particular the sign of Ce is not determined by perturbative cosidera- 
tions alone. The latter requires additional non-perturbative knowledge of matrix elements 
of composite operators. A crude estimation using the non-relativistic quark model indi- 
cates that Ce is positive, which implies a repulsive core with a potential diverging a little 
weaker than the generically expected r~ 2 at short distances. 



2. Operator Product Expansion and potentials at short distance in 3 di- 
mensions 

In this section we discuss the application of the operator product expansion (OPE) to the 
determination of the short distance behavior of the BS potential. We consider the equal 
time Bethe-Salpeter (BS) wave function defined by 

<p%{r) = (0\O A (r/2, 0)O B (-f/2,0)\E) , (2.1) 

where \E) is an eigen-state of the system with energy E, and A , Ob are some operators 
of the system. Here we suppress other quantum numbers of the state \E) for simplicity. 
Using the OPE of A and O b 

O A (f/2,0)O B (-f/2,0) ~ J2 D Ab(^)O c (0,0) , (2.2) 

c 

1 Unfortunately the results for /3 01 and f3 10 as given in ref. & differ (incorrectly) from ( [L.4| ) by a factor 
of 2. 
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we have 

^s(r) * ^^ B (r)(0|Oc(0,0)|E) . (2.3) 

c 

Note that r dependence appears solely in D^ B (f) while the E dependence exists only in 
(0\Oc(0, 0)\E). As we will see, in the r = \r\ — > limit, the coefficient function behaves as 

DabW - r ac (-\ogrfcfc{eA) , (2-4) 
where 9, <j) are angles in the polar coordinates of r, so that 

PAB(f) ^X>*(-logr)^/ c (M)I?c(£), I>c(S) = (0\O c (0,0)\E) . (2.5) 

We now assume that C has the largest contribution at small r: 

«c < ac or (2.6) 
a c = a c >, Pc > Pc ■ (2-7) 

for V C" / C. The potential can be calculated from this wave function. 

As will be seen later, ac = o-c 1 = for the NN case in QCD. Furthermore states with 
zero orbital angular momentum [L = 0) dominates in the OPE, so that the wave function 
at short distance is given by 

<PAu(r) - [(-\ogr)fcD c (E) + (-logr)^'£> c ,(£)] (2.8) 

with /3c > Pa- Using 

V 2 (-logr)^ = -p{-\ogrf~ l r~ 2 , (2.9) 

— log r 

we obtain the following classification of the short distance behavior of the potential. 



1. Pc / 0: The potential at short distance is energy independent and becomes 

V(r) ~ - ^ (2.10) 
r z (— logr) 

which is attractive for Pc > and repulsive for Pc < 0. 

2. Pc = 0: In this case we have 

The sign of the potential at short distance depends on —Pc'Dc(E)/Dc(E). 

If there are two or more operators which have the largest contribution at short distance, 
we have 

VabW = (-logr)te(D Cl (E)+D C2 (E) + ■■■)■ (2-12) 

In this case, the above analysis can be applied just by replacing Dc{E) — > Dc 1 (E) + 
D C2 (E) + ---. 

On the lattice, we do not expect divergences at r = due to lattice artifacts at short 
distance. The above classification holds at a <C r < 1 fm, while the potential becomes 
finite even at r = on the lattice. 
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3. Renormalization group analysis and operator product expansion 

3.1 Renormalization group equation for composite operators 

In QCD, using dimensional regularization in D = 4 — 2e dimensions, bare local composite 
operators (x) are renormalized [@] according to 2 

{ p\x) = Z AB {g,e)Of{x). (3.1) 

(Summation of repeated indices is assumed throughout this paper.) The meaning of the 
above formula is that we obtain finite results if we insert the right hand side into any 
correlation function, provided we also renormalize the bare QCD coupling go and the 
quark and gluon fields. For example, in the case of an n-quark correlation function with 
operator insertion, which we denote by G^ A {go,e) (suppressing the dependence on the 
quark momenta and other quantum numbers) we have 

g£lW) = Z AB (g,e)Z- n/2 (g,e)g%(g ,e). (3.2) 

We recall from renormalization theory that for the analogous n-quark correlation function 
(without any operator insertion) we have 

Gt Q) (9^) = Z- n/2 (g,e)g^(g ,e), (3.3) 

where the coupling renormalization is given by 

g 2 =^Z 1 (g,e)g 2 . (3.4) 

The renormalization constant Z\ in the minimal subtraction (MS) scheme we are using has 
pure pole terms only: 



Z 1 {g,e) = l — + -y_+O(0°), 3.5) 

2e e z 



where 



16vr 2 \ 3 3 s ) 256vr 4 \ 3 \3 N J 1 J y J 

Similarly for the fermion field renormalization constant, we have 

Z F (g,e) = l-^ + 0(g% (3.7) 

where jfq is given by ( 4.1C| ). The gluon field renormalization constant is also similar, but 
we do not need it here. Finally the matrix of operator renormalization constants is of the 
form 

(1) 2 

Z AB (9,e) = Sab - + 0( 5 4 ). (3.8) 



2 We note that we are considering the massless theory here since quark masses play no role in our analysis. 



- 5 - 



The renormalization group (RG) expresses the simple fact that bare quantities are 
independent of the renormalization scale [i. Introducing the RG differential operator 

the RG equation for n-quark correlation functions can be written as 

{v + ^ F (9)}Gt n) (9^)=0, (3.10) 
where the RG beta function is 

P{g) = eg + (3 D (g, e) = eg - - = -fog 3 - Pl9 5 + 0(g 7 ), (3.11) 

1 + 2 8g 

where /3r>(g, e) is the beta function in D dimenions and the RG gamma function (for quark 
fields) is 

Ma) = p D (g,e) d -^- = jFog 2 + 0( 5 4 ). (3.12) 

It is useful to introduce the RG invariant lambda-parameter A by taking the Ansatz 

A = /;e /(9) (3.13) 

and requiring DA = 0. The solution is the lambda-parameter in the MS scheme (Ams) if 
the arbitrary integration constant is fixed by requiring that for small coupling 

/(9) =-^-4 l " ( "° 92)+0(£ ' 2) ' (3 - i4) 

Finally the RG equations for n-quark correlation functions with operator insertion are of 
the form 

{v + Iif{ 9 )} G^ ) (g,v)-jAB(g)G { :™\g^)=0, (3.15) 

where ^ 

lAB {g) = -Z AC (3 D (g, e)^|p = lA %g 2 + 0(/). (3.16) 

3.2 OPE and RG equations 

Let us recall the operator product expansion ( [2.2j ) 

0!(y/2)0 2 (-y/2) ~ D B (y) O fl (0). (3.17) 

We will need it in the special case where the operators 0\ , O2 on the left hand side are 
nucleon operators and the set of operators Ob on the right hand side are local 6-quark 
operators of engineering dimension 9 and higher. All operators in (3.17) are renormalized 
ones, but from now on we suppress the labels ( ren ). As we will see, the nucleon operators 
are renormalized diagonally as 

1 =( 1 (g,e)0{° ) , 2 = C, 2 {9,z)0f\ (3.18) 
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and we can define the corresponding RG gamma functions by 

71,2(5) = Pd( 9 , e) = 7i> 2 + 0(g 4 ). (3.19) 

Next we write down the bare version of fl3.17|) (in terms of bare operators and bare 
coefficient functions): 

o[ 0) (y/2)O { °\-y/2) ~ D<g\y)O$\0). (3.20) 

Comparing ( gl7| ) to ( jDoD , we can read off the renormalization of the coefficient functions: 

D B (v) = Ci(9,e)C 2 (9,e)Df(y)Z A 1 B (g,e) (3.21) 

and using the /U-independence of the bare coefficient functions we can derive the RG equa- 
tions satisfied by the renormalized ones: 

VD B (g, n, y) + D A (g, fx, y) j AB (g) = 0, (3.22) 

where the effective gamma function matrix is defined as 

1ab{q) = iAs{g) - [71(5) +72(5)] S A b- (3.23) 

3.3 Perturbative solution of the RG equation and factorization of OPE 

We want to solve the vector partial differential equation (|3.22j ) and for this purpose it is 
useful to introduce Uab{9), the solution of the matrix ordinary differential equation 

P{g) U AB (g) = iAc{g) U CB (g) (3.24) 
dg 

and its matrix inverse Uab (fiO- We will assume that the coefficient functions are dimen- 
sionless and have the perturbative expansion 

D A (g, /x, y) = D A (g; M r) = D A . + g 2 D A]1 ^r) + 0( 5 4 ), (3.25) 

where r = \y\. For the case of operators with higher engineering dimension 9 + a the 
coefficients are of the form r a times dimensionless functions and the analysis is completely 
analogous and can be done independently, since in the massless theory operators of different 
dimension do not mix. (In the full theory quark mass terms are also present, but they 
correspond to higher powers in r and therefore can be neglected.) 

We will also assume that the basis of operators has been chosen such that the 1-loop 
mixing matrix is diagonal: 

7 ab(s) = 2^(3 A g 2 5 AB + 0(g i ). (3.26) 



In such a basis the solution of ( 3.24 ) in perturbation theory takes the form 

UAs(g) = {Sab + R AB (g)} <T 2fe , (3.27) 
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where Rab{9) = 0(g 2 ), with possible multiplicative log g 1 factors, depending on the details 
of the spectrum of 1-loop eigenvalues /3a- 



Having solved (3.24) we can write down the most general solution of ( |3.22j ): 

D B (g;fir)=F A (Ar)UAB(g)- 



(3.28) 



Here the vector Fa is RG-invariant. Introducing the running coupling g as the solution of 
the equation 

f{g) = fig) + MH = h(Ar) 



Fb can be rewritten as 

F B {Ar) = D A {g-l)U AB {g). 
Since, because of asymptotic freedom (AF), for r — > also g — > as 



2/?oln(Ar)' 



(3.29) 
(3.30) 

(3.31) 



Fb can be calculated perturbatively using ( |3.25| ) and ( p. 27 ). 

Putting everything together, we find that the right hand side of the operator product 
expansion ( 3.17] ) can be rewritten: 



where 



Oi(y/2)O a (-y/2) ~ F B (Ar)O B (0), 



B = U BC (g)O c . 



(3.32) 



(3.33) 



There is a factorization of the operator product into perturbative and non-perturbative 
quantities: F B (Ar) is perturbative and calculable (for r — > 0) thanks to AF, whereas the 
matrix elements of B are non-perturbative (but r-independent). 

An operator O b first occurring at ^e-loop order on the right hand side of ( |3.17| ) and 
corresponding to normalized 1-loop eigenvalue (3 B has coefficient F B (Ar) with leading short 
distance behavior 



F B (Ar) 



9 



2((b-/3b) 



-2A>ln(Ar)) 



(3.34) 



In principle, an operator with very large f3 B , even if it is not present in the expansion 
at tree level yet, might be important at short distances. This is why it is necessary to 
calculate the full 1-loop spectrum of f3 B eigenvalues. As we shall see, no such operators 
exist in our cases, and therefore operators with non-vanishing tree level coefficients are 
dominating at short distances. The corresponding coefficient functions have leading short 
distance behavior given by 



F B (Ar)wD B ;o (-2/? ln(Ar)) 



0B 



(3.35) 



A similar analysis in the case of operators of dimension 9 + a leads to the result (| 
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4. OPE and Anomalous dimensions for two nucleons 



4.1 OPE of two nucleon operators at tree level 

The general form of a gauge invariant local 3-quark operator is given by 

Bf(x) ee B^x) = e« bc q^{x)qf{x)qf{x) , (4.1) 

where a, (3, 7 are spinor, /, g, h are flavor, a, b, c are color indices of the (renormalized) quark 
field q. The color index runs from 1 to N = 3, the spinor index from 1 to 4, and the flavor 
index from 1 to Nf. In this paper a summation over a repeated index is assumed, unless 
otherwise stated. Note that Bj 1 * is symmetric under any interchange of pairs of indices 
(e.g. = BpJ^) because the quark fields anticommute. For simplicity we sometimes 

use the notation such as F = fgh and F = afi^ as indicated in flO]). 
The nucleon operator is constructed from the above operators as 

Bf(x) = (P + ,) aa , B f J h pi {C l5 ) M (inY h , (4-2) 

where P+4 = (1 + 74)/2 is the projection to the large spinor component, C = 7274 is the 
charge conjugation matrix, and T2 is the Pauli matrix in the flavor space (for Nf = 2) 
given by (iT 2 ) f9 = e /9 . Both C75 and it 2 are anti-symmetric under the interchange of two 
indices, so that the nucleon operator has spin 1/2 and isospin 1/2. Although the explicit 
form of the 7 matrices is unnecessary in principle, we find it convenient to use a (chiral) 
convention given by 



Ik = . „ , 74 = L, n , 75 = 71727374 = n . • (4.3) 






As discussed in the previous section, the OPE at the tree level (generically) dominates 
at short distance. The OPE of two nucleon operators given above at tree level becomes 

Bf a (x + y/2)B9(x - y/2) = Bfa)B«(x) + ^ [d^B^B^x) - B^d^x)]} 



+ 



[dMBl(x)B 9 p (x)] - 4d„Bl(x)d„B°(x)} + ■■■ .(4.4) 



For the two- nucleon operator with either the combination [a/3], {fg} (S = 0) or the com- 
bination {a/3}, [fg] (S = 1), terms odd in y vanish in the above OPE, so that only 
even L contributions appear. These 6-quark operators are anti-symmetric under the ex- 
change (a,/) (f3,g). On the other hand, for two other operators with ([a/3], [fg]) or 
({a/3}, {fg}), which are symmetric under the exchange, terms even in y vanish in the OPE 
and only odd L's contribute. 

Knowing the anomalous dimensions of the 6-quark operators appearing in the OPE, 
which will be calculated later in this section, the OPE at short distance (r = \y\ <C 1, 
y4 = 0) becomes 

B{{x + i//2)BJ(x - y/2) ~ £ c A (r)O f a f(x) + £ d B (r)y k y l O f a f kl (x) 

A B 

+ Ys e c(r)y k O{f k (x) + --- , (4.5) 

c 
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where the coefficient functions behave as 

ca(t) ~ (— logr)^ , dsir) ~ (— logr)^ 3 , ec(r) ~ (-logr)^ , (4.6) 

and Pa,b,c are related to the anomalous dimensions of the 6-quark operators Oj^ ,A , of 
those with two derivatives O and of those with one derivative 0^ k . 
The wave function defined through the eigenstate \E) is given by 

vT n (v) = (0\Bl(x + y/2)Bl(x-y/2)\E) 

^ -£c A (r)(0\Olf(x)\E) +Y / dB(r)y k y l {0\O f c $ B kl (x)\E) + 0(/) , (4.7) 

A B 

for the anti-symmetric states, while 

C 



<p° E dd (y) = (0\Bf a (x + y/2)B°(x - y/2)\E) ~ £ e c (r)y fc (0|O^(x)|E) + 0(y 3 ) (4.8) 



for the symmetric states. 

In this paper, we consider only 6-quark operators without derivatives and calculate 
the corresponding anomalous dimensions. 

4.2 General formula for the divergent part at 1-loop 

Following the previous section, we define the renormalization factor Zx of a /c-quark op- 
erator X = [q k ] through the relation 

[q k ]™ = Z x [q k ] = Z x Z k J 2 [q k ], (4.9) 

where qo(q) is the bare (renormalized) quark field. The wave function renormalization 
factor for the quark field is given at 1-loop by 

z F = i + ?zU, z? = --^r e (4.io) 

where A is the gauge parameter and C F = N ^ x l • 

At 1-loop the renormalization of simple fc-quark operators (those involving no gauge 
fields) is given by the divergent parts of diagrams involving exchange of a gluon between 
any pair of quark fields. The 1-loop correction to the insertion of an operator q& ^ (x)q b ^ 9 (x) 
in any correlation function involving external quarks is expressed as the contraction of 

Qa J (x)q b f(x)^ j d D yd D z^(y)^(z)[^Ky)^T A 7^ / Hy)][^ 1 ^)^T i? 7^ 9l (^)](4.11) 

where tr T A T B = 5 A b/2 

in our normalization. Since two identical contributions cancel the 
2! in the denominator, the contraction at 1-loop is given by 

-g 2 {T A ) aai {T A ) bbl j d D yd D z [S F (x - y)^q{y)] a a ' h G, u (y - z) 

x [S F (x - zh u q(z)] b fs 191 (4.12) 



- 10 - 



where the free quark and gauge propagators are given in momentum space as 



S F (p) 



ip' + m 



p 2 + m 2 



A: 2 



V - (! - A ) 



k^k u 



(4.13) 



The above contribution can be written as 

,2 r d D pd D q 

|2£> 



(27T 



Taa 1>m (P, q) C/ 1 (P)e^gjf (g)e^ (4.14) 



where 



T aai ^ h (p,q) = J -j—^ [SF{p + k)~iJ\ aai G^ u {k) [S F {q - k)^]^ 
whose divergent part is independent of the momenta p, q and is given by 



(4.15) 



?aai,A8i(0, 0) 



1 1 



16vr 2 e 



4£. 



<V/ (8 cx Mi , + Al <8> 1 



(4.16) 



with cr^ = | [7^,7^]- We then obtain the divergent part of the 1-loop contribution as 



q a /(x)q b f(x) 



1— loop,div 



2 1 f 

''(To + AT, ) • q a (x) q\x 



32iVvr 2 e 

where (bold-faced symbols represent matrices in flavor and spinor space) 



a,/3 



(Ti)^ 1 ^ = {1 ® 1 + iVKgl}^ 1 ' 9 ^ . 
Here we use the notation 

{X Y}£^ = X& Yg£ {XSY}^ = Xj* , 

Using the following Fierz identities for spinor indices 

~7 £ °> ® ^ = Pr® Pr + Pl® Pl- 2{P r ®P r + P l ®Pl) , 



7 £ <V®<V = Pr®Pr + Pl®Pl - 2(Pr ®Pr + P l ® Pl) , 



where Pr, Pl are the chiral projectors 



Pr = ^(1 + 75), Pl = ^(1-75) 



(4.17) 

(4.18) 
(4.19) 

(4.20) 
(4.21) 

(4.22) 
(4.23) 

(4.24) 



we can simplify To as 

( T o) f JaM = 6 " 16991 [<Wa8i " 2 <W<Wi] + N5^5^ [5p ai 5 aPl - 25 aai 



(4.25) 



where either ot\,fii £ {1, 2} (right-handed) or a\^\ 6 {3, 4} (left-handed) due to the chiral 
projections in eqs. ( 4.22j ) and ( [4.231) . In our following calculation of the 1-loop anomalous 
dimensions, eq. ( |4.17 ) together with eqs. ( |4.25 ) and (4.19) are the key equations. 
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4.3 Renormalization of local 3 quark operators at 1-loop 

In this subsection we calculate the anomalous dimensions of general 3-quark operators at 
1-loop. In terms of the renormalization factor defined as 

C[q 3 ] = (z 3 J 2 [ q % c = i + g 2 (C (1) + Cf) + • • • , 



Tjrcnor. 



(4.26) 



where ( ) 18 the A-independent (dependent) part at 1-loop, the divergent part of 
the insertion of the 3-quark operator Bf = B*^f* defined in ( |4.1| ) at 1-loop is given by a 
linear combination of insertion of baryon operators, and (with a slight abuse of notation) 
we express this as 



( r (i)div)F 



z / rr' 



The A-dependent contribution from Ti in ( 4.19 ) is diagonal and given by 

9 2 N + l F 



2/ r (l)divsF 



3A- 



-B 



(4.27) 



(4.28) 



so that the A-dependent part of Q vanishes: 



Ai) 



3A N + l 3 (1) 
32vr 2 Ne 2 



A 3(iV + l)(JV-3) 



0. 



(N = 3). (4.29) 



64/Vtt 2 e 

Therefore ( is A-independent, as expected from the gauge invariance. We remark that we 
leave A?" explicit in some formulae to help keep track of the origin of the various terms, but 
in our case we should always set = 3 at the end. 

The A-independent part of I^ 1 ) from To in ( 4.25 ) leads to (N = 3): 

(N + l) g 2 



/■p(l)div\fgh 

f r (l)dw\fgh 
v 1 taffi 



orqfgh 

2N 167r 2 e L a ^< 
{N + l) g 2 



9 f>f9h _ 9 rjfgh _ 9 rjfgh 
Z " D /3q7 LD ~iPa LD a~jP 



gfgh _ 2B^ gh 

2N 167r 2 e L a ^ 



(4.30) 
(4.31) 



where a,/3,7 G {1,2} (right-handed), while 7 G {1 = 3,2 = 4} (left-handed). Note that 
the same results hold with hatted and unhatted indices exchanged. 

These relations can be easily diagonalized and the combinations which do not mix are 
given by 



(AV)f9h 
^ ){aa/3} 


- (A~t)\}9 h - 
vs 


€ 


(4.32) 


(Athf^gh 

^ )[aP]a 


_ /Aihfrgh - 
vs '[&$]& 


d 

-12- , 


(4.33) 


(AW 9h 


_ /Aihfgh _ 




(4.34) 


(Al))&9h 


_ (Ai))f^gh _ 


e 


(4.35) 



where d is given by 



1 



1 



32Ntt 2 96vr 2 



(4.36) 
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The square bracket denotes antisymmetrization [a/3] = a/3 — f3a, and curly bracket means 
{a/3} = a/3 + f3a, {aa/3} = aa/3 + a/3a + j3aa. The totally symmetric case corresponds to 
the decuplet representation (for Nf = 3) and contains the Nf = 2,1 = 3/2 representation. 
The antisymmetric case corresponds to the octet representation (for Nf = 3) and contains 
the Nf = 2 , 1 = 1/2 representation. 

The anomalous dimension at 1-loop is easily obtained from 

7 = 5 2 7 (1) + 0( 5 4 ) = (3 D (g, e)^ = -2(^g 2 e + 0( 5 4 ) . (4.37) 



Therefore we have 



( 7 d; 



} = 7 ( n - = -24d, 4.38 

/ {aa/3} V / {&a/3} 

U^ 9h =U^ 9h = 2U, (4.39) 

V / [a/3] a V / [60\& 

^)) f9h =U^) f9h =-8d, (4.40) 
Ui)) M9h =Ui)) M9h = 24d. (4.41) 

V / [a/3] 7 V / [a/3] 7 

4.4 Anomalous dimensions of 6 quark operators at 1-loop 

In this subsection we consider the renormalization of arbitrary local gauge invariant 6— 
quark operator of (lowest) dimension 9. Any such operator can be written as a linear 
combination of operators 

O c {x) = B^f 2 (x) = B^(x)Bg(x) = O a {x)O b {x) , (4.42) 

with A = (Tf, Ff) and B = (T2, F2). Note Ojf[x) and/or Ob(x) may not be operators with 
proton or nucleon quantum numbers and separately may not be diagonally renormalizable 
at one loop. The reason for considering the renormalization in more generality is that 
in principle there may be operators in this class which occur in the OPE of two nucleon 
operators at higher order in PT, but are relevant in the analysis because of their potentially 
large anomalous dimensions. 

4.4.1 Linear relations between 6 quark operators 

According to the considerations in subsect. 4.2 the operators in eq. Q4.42D mix only with 
operators Oc = Oa'Ob' which preserve the set of flavors and Dirac indices in the chiral 
basis i.e. 

FfUF 2 = f{ u f' 2 , Tf u r 2 = r' x u . 

Note however that such operators are not all linearly independent. Relations between them 
follow from a general identity satisfied by the totally antisymmetric epsilon symbol which 
for N labels reads 

jY £ a i--- a JV £ bi...b N _ ~y ' £ a 1 ...aj^ 1 b k aj +1 ...a N £ b 1 ...b k _ 1 ajb k+1 ...b N (4 43) 
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For our special case, N = 3, this identity implies the following identities among the 6-quark 
operators 

3Br F !;r F ;+E<;Sii = 0. (4-44) 

where i-th index of abc and j-th index of def are interchanged in (abc,def)[i,j]. For 
example, (ri,r 2 )[l,l] = a 2y 5i7i, ai/?272 or (ri,r 2 )[2, 1] = ai« 2 7i, /3i/3 2 7 2 . Note that the 
interchange of indices occurs simultaneously for both T\ , T 2 and F\ , F<i in the above formula. 
The plus sign in ( 4.44j ) appears because the quark fields are Grassmann. 



An immediate consequence of the identity is that the divergent part of the A-dependent 
contributions, calculated from Ti in ( 4.17| ), must vanish, after the summation over the 9 



different contributions from quark pairs on the different baryonic parts A, B is taken. The 
A-dependent part of the contribution of quark contractions on the same baryonic parts 
is compensated by the quark field renormalization. Thus the renormalization of the bare 
6-quark operator is A-independent as expected from gauge invariance. 

As an example of identities, we consider the case that Ti,r 2 = aa(3,af3f3 (a ^ (3 and 
Fi,F 2 = ffg, ffg (/ ^g), the constraint gives 

od/ZsJ/j i (o _ <)\ -nffgjfg , p>fff,fgg , /n _ -i \ -nfggjff 
_ Ar>ffg,ffg 1 rifffjgg , nfggjff _ n (a ak\ 

where minus signs in the first line come from the property that B^'rl = ~-^rl'rt' There 



are no further relations among 6-quark operators beyond ( 4.44 ). 
4.4.2 Divergent parts at 1-loop 

We thus need only to calculate the contributions from To, which can be classified into the 
following 4 different combinations for a pair of two indices: 

A (A -*-(N +!)(()({) , (4.46) 



j) =► (1 - 2N) (/) + (N - 2) (>) (/) , (4.47) 

£)(£) ^ -(£)(£) -*(£)(£) , (4.48) 

(£) (?) - {(*) (?) - 2 (?) (£)}+"{(?) (£) - 2 (£) (?)} . ( 4 - 49 ) 

where f ^ g and a ^ j3 E (1,2) (Right) or E (3,4) (Left). 

The computation can be made according to the following steps: 

i. ) Select the total flavor content e.g. 3/ + 2>g or 4/ + 2g (f ^ g). These are the only 
cases we will consider since in this paper we are mainly restricting attention to baryon 
operators with TVf = 2, but the approach is also applicable to more general cases (iVf > 2). 

ii. ) Given a flavor content classify all the possible sets of Dirac labels in the chiral 
basis e.g. 111223,112234,... It is obvious from the rules above that some have equivalent 
renormalization at 1-loop e.g. 111223 and 112223 with 1 -H- 2, and also those with hatted 
and unhatted indices exchanged e.g. 111223 and 133344. 
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iii. ) For given flavor and Dirac sets generate all possible operators 3 . Then generate all 
gauge identities between them and determine a maximally independent (basis) set {Si}. 

iv. ) Compute the divergent parts of the members of the independent basis: 



v.) Finally compute the eigenvalues and corresponding eigenvectors of ^ T to determine 
the operators which renormalize diagonally at 1-loop. 

An example of the procedure is given in Appendix A. Some of the steps are quite 
tedious if carried out by hand. e.g. in the case 3/ + 3g and Dirac indices 112234 there 
are initially 68 operators in step iii. with 38 independent gauge identities, and hence an 
independent basis of 30 operators. However all the steps above can be easily implemented 
in an algebraic computer program using MATHEMATICA or MAPLE. 

If the quarks /, g belong to an iso-doublet e.g. we identify / with u having /3 = 1/2 and 
g with d (I3 = —1/3), then if an eigenvalue is non-degenerate the corresponding eigenvector 
belongs to a certain representation of the isospin group. If the eigenvalue is degenerate 
then linear combinations of them belong to definite representations. For the 3/ + 3g case 
they can have I = 0,1,2,3. Eigenvectors with / = 0,2 are odd under the interchange 
/ -H- g and those with 7 = 1,3 are even. The operators in the case 4/ + 2g have 13 = 1 
and hence have I = 1,2, 3. The eigenvectors in this case can be obtained from those of the 
3/ + 3g case by applying the isospin raising operator. 

The complete list of eigenvalues and possible isospins are given in Tables Q4l| in Ap- 
pendix A. Here we summarize the most important results. 

1) For the 3/ + 3g (and 4/ + 2g) cases all eigenvalues jj < 48d = 27^, where 77V is 
the 1-loop anomalous dimension of the nucleon (3-quark) operator. We have not found an 
elegant way of proving this other than computing all cases explicitly. 

2) It is easy to construct eigenvectors with eigenvalue 277V e.g. operators of the form 
B,a i-B - - since there is no contribution from diagrams where the gluon line ioins quarks 
in the different baryonic parts. 

3) Operators with higher isospin generally have smaller eigenvalues. 

4.4.3 Decomposition of two-nucleon operators 



idiv 
i 




(4.50) 



Since 



C75 



/o -1 \ 
10 



(4.51) 



-1 



\0 1 / 



in the chiral representation, the nucleon operator is written as 



a+a,[2,l] 




(4.52) 



3 recall the single baryon operators are symmetric under exchange of pairs of indices 
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where a = 1,2, a = a + 2, and / 7^ g. This has anomalous dimension 7jv = 2Ad. 

We then consider two independent 6-quark operators occurring in the OPE at tree 
level which decomposed as follows. The spin-singlet (S = 0) and isospin-triplet (/ = 1) 
operator is decomposed as 



jffg gffg 

a+a,[p,a]+\p,&] /3+$,[P,a]+0,a] 



B J ' " x B{ JS : ,„ , rS = Bf + Bfj + Bf H + B% • B° ] • />><", 



(4.53) 



where 



B 01 



B oi 



1 

01 
11 



B 



01 



B 



01 



R oi 

R 01 



//9 



a[/3,a: 
ce[/3,a 

a[/3,ce 

= B^? 

a[/3,a 



ffg 

a[/3,a 



B ff 9 

gffg 
gffg 



B 



ffg 
ffg 



B 



+ B 
+ B 
+ B 
+ B 
+ B 
+ B 



ffg r>ffg 



ffg 

a\f},& 

ffg 

a[/3,a 

ffg 

d[/9,a 

ffg 

d[/9,a 

ffg 

&[$,& 



B 

$[$,&] ' 

r//9 1 r//# r//9 1 r//S r//9 



Tiff 9 i_ r//s 

^/3[/3,a] + ^ 



,1? 



a[/3,a] /3[/3,« 



+ 5 



//9 



Tiff 9 



B 



ffg 



B 



ffg 



(4.54) 
(4.55) 
(4.56) 
(4.57) 
(4.58) 
(4.59) 



where a / In the above we do not have to calculate all contributions. Some of them 
are obtained from interchanges under (1,2) o (3,4) or (1,3) o (2,4). 

Similarly the spin-triplet (S = 1) and isospin-singlet (I = 0) operator is decomposed 

as 



gffg ^ggf 

a+a,[0,a] + [l3,a} a+a,\fi,a]+[j3,a] 



r10 1 r10 1 r10 1 R 10 1 r10 1 r10 
B I + II + III + B IV + By + By j , 



(4.60) 



where 



10 



B 



r10 
B II 

r10 
III 

rIO 



13 



10 



RlO 



B 



ffg gggf 



B 



= B 



a[/3,a 

ffg 

a[l3,a 

ffg 



fiffg 

a[/3,a 

gffg 



B 



ffg 



B 



a[/3,a 

ggf 
ggf 



ffg gggf 



B 

"LP, 

B"t 

a[/3,& 

B 99 ,{ 

a[p,a 



B 



ggf 



a\p,a 

+ B ff ? . 

aSfi,6t 
+ &\p,a 



+ B 



ffg 

a[/3,a 



gffg 

a[/3,a 



+ B 



ffg 

&\p,& 



&[$,& 
ggf 

a\p,a 

ggf 

' a\p,a 



B 



B 



B 9 Vr 
a[p,a 

B"t 



B 



ggf 



+ B ff ? B 99 { , + B f J 9 iB 99f 



_l_ gffg b"1 

&[$,&] a[0,a) 



afoot] afoot] 

_l_ £>ffg B"i 

otfoa] otfoa] 



(4.61) 
(4.62) 
(4.63) 
(4.64) 
(4.65) 
(4.66) 



Again a half of the above 1-loop contributions can be obtained from others by the inter- 
change (1,2) o (3,4) or /-> g. 



AAA Results for anomalous dimensions 

It is very important to note here that operators By^ for both cases (SI = 01 and 10) have 
the maximal anomalous dimension at 1-loop, since as noted in point 2) above, no 1-loop 



correction from Tq joining quarks from the two baryonic components exists for B 



F U F 2 
af3f,a' $'7' 
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type of operators. Therefore we always have some operators with (3a = which dominate 
in the OPE at short distance. 

The 1-loop corrections rW to 6-quark operators B SI are computed in appendix [A| and 
are summarized as: 

(rf) W = -12~Bf , (r?}) (1) = 12^ 7 °) , (Tfjj) (1) =0, 

(r%) {1) = 0, (TVf } =6^ + 6^, (r«) (1) =24^, (4.67) 

for SI = 01. The last two results can be written as 

(I$) W = 6-B° v \ , (r y V) (1) = 24^<, , , (4.68) 



where 



B$ = By 1 - , B v \, = fl& . (4.69) 



Similarly we have for SI = 10 



(r?) (1) = -4V, (r)?) (1) = 2oi$, (r]? 7 ) (1) = o, 

(rf v ) (1) = 8±B% , (rt?,) (1) = 6^B# , = 24^#» , (4.70) 



where 



r10 _ r10 r10 r 10 _ rIO o 7-1 \ 

By, — By - ^Byj , By j, — Byj . (4. t L) 

Denoting the eigenvalues of the anomalous dimension matrix by 7c, we give the values 
of 7 s7 defined by 

7C - 2 1N = 2dj SI , (4.72) 

in table [l] (N = 3), which shows, in both cases, that the largest value is zero while others 
are all negative. The case 2 in sect. |2| is realized: /3c = and 

/3c'=/3o 01 = - 33 _ 6 2jVf for5 = 0,/ = l, (4.73) 

^=/?o 10 = - 3 3_ 2 2Arf for 5 = 1,1 = 0. (4.74) 



5. Short distance behavior of the nucleon potential 

We consider the following structure of the potential. 

V(y) = V (r) + F CT (r)(<ri • a 2 ) + V T (r)S 12 + O(V) (5.1) 

where r = \y\, and 

S\ 2 = ^{a x ■ y)(a 2 ■ y) - [a\ ■ a 2 ), V = tL (5.2) 

\y\ 
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Table 1: The value of 7 s7 (denned in (4.72)) for each eigen operator in the SI = 01 and SI = 10 
states. 





1 


77 


m 


JV 


V 


VI' 




-36 


-12 


-24 


-24 


-18 







-28 


-4 


-24 


-16 


-18 






is the tensor operator. Here d{ acts on the spin labels of the i th nucleon. 

Since, as shown in the previous section, 6-quark operators appeared at tree level in 
the OPE of NN which have the largest and the second largest anomalous dimensions, we 
mainly consider 6-quark operators at tree level in the OPE, which is written as 

B{{x + y/2)B*(x - y/2) ~ cvjB^^x) + c n (- \ogrf"' B{ 9 Ia ^x) + ■■■ (5.3) 

where cyi and c/j are some constants, and • • • represents other contributions, which are 
less singular than the first two at short distance. (We here write the spinor and flavor 
indices a, ft and /, g explicitly for later use.) The anomalous dimensions ft^ 1 are given in 
( 14731) and ( gj7g ). 

5.1 Potential for 5 = and 1=1 states 

In the case that S = and I = 1, we take a ^ ft and / = g in eq. (]5.3| ) , whose leading 
contributions couple only to the J = L = state, which is given by 

\E) = \L Z = 0,S Z = 0, I z = l) L =o,s=o,l=i = |0, 0, l) 0)0) i • (5.4) 

The relevant matrix elements are given by 

^1(01^^10,0, l)o,o,o = AyjY^aftKfgh, (5.5) 



c//<0|s// iaj9 |0,0, l)o,o,o = A Q n Y^[aft]{fg} u (5.6) 



where A^j and Ayj are non-perturbative constants, Y^ z is a spherical harmonic function, 
[aft] = (5 a i5p2 — SpiS a 2)/V2 represents the (S,S Z ) = (0,0) component, and {fg}i = <5/i^gi 
corresponds to isospin (I,I Z ) = (1,1). With the notation that <fi( 1 So,J z = 0) IIz=11 = 
YQ[aft]{fg}i, the wave function at short distance is dominated by 

^ E So (y) = (0\Bf(z + y/2)B°(x- y/2)|0, 0, l) 0)0 ,i 

' Ayj + A° H {- logr)^ 01 ) So, 0) 11 + • • • (5.7) 



from which we obtain 

^^ ( -'° e f"-^- ^0)" + ... (5.8) 
2m r z mjv 

where m = mjv/2 is the reduced mass of the two nucleon system. 
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Since S12 is zero on 0( 1 S'o,O) 11 , we have 

K 01 (r)^°(y) 



where V^ 01 (r) = Vo(r) — 3V r f7 (r). We therefore obtain 



V^{r) ~ F Ui (r) 



01, 



/7 



4° ' 



where 



F s/ (r) 



(5.9) 



(5.10) 



(5.11) 



The potential diverges as F (r) in the r — > limit, which is a little weaker than r . 
5.2 Potential for 5 = 1 and 1 = 

We here consider the spin-triplet and isospin-singlet state(5' = 1 and 1 = 0). Since 1 = 
and I z = in this case, we drop indices I and I z unless necessary. In this case, the leading 
contributions in eq.(|5.3|) couple only to the J = 1 state, which is given by 



where 



\E) = | 3 Si, J z = l)+x\ 3 D 1 ,J z = 1), 



3 Si, Jz = 1) = \L Z = 0, S z = l)i=o,s= 



5 £>i,J z = l) 



|0,1)-V3|1,0) + a/612, -1) 



L=2,5=l 



(5.12) 

(5.13) 
(5.14) 



x is the mixing coefficient, which is determined by QCD dynamics, and 2S+1 Lj specifies 
quantum numbers of the state. 

Relevant matrix elements are given by 



Ci(0 



D 



fy 



5i,l) =B°4>( 3 S 1 ) , 



Ci(0 



0. 



for i = II and V/, where are non-perturbative constants, and 

<f>^S x ) =y o °(0,0){a/3} 1 [/ 5 ], 

with {a/3}i = Jai^i. 

Using the above results, we have 

v J E =1 (y) - {^/ + (-logr)^^ / }^(3 5l ) 



(5.15) 



(5.16) 



By applying V 2 we obtain 



V 



Si) 



(5.17) 



(5.18) 
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while 

V(y)ip J E =1 (y) ~ B^V^ir)^ ( 3 5i) + 2V2Fr(r)B£ J ( 3 Di) , (5.19) 
where we use the formula in the appendix [B|, V^ 10 (r) = Vq(t) + V^(r), and 

{ 3 Di) = -L [y 2 °{a/3}i - V^H^Io + \/6Y 2 2 {a/3}-il [/</] , (5-20) 



with {a/3}_i = 5q 2 5/32 and {a/3} = (S a i$/32 + 5^6^)/^. 
By comparing eq. Q5.18p with eq.( |5.19| ), we have 

B° TT 



V c 10 (r) ~ F 10 (r)-f , F T (r) ~ . (5.21) 
% 

This shows that the central potential V^, 10 (r) diverges as i ?10 (r) in the r — > limit, which 
is a little weaker than 1/r 2 , while the tensor potential Vr(r) becomes zero in this limit at 
the tree level in the OPE. 

5.3 Higher order in the OPE and the tensor potential 

While the short distance behavior of the central potential is determined by the OPE at the 
tree level, the determination of the tensor potential at short distance requires the OPE at 
higher order, whose relevant contribution is given by 

Bi(x + y/2)B 9 (x - y/2) ~ c VI B f v 9 I ap {x) + c n (- \ogrf"' B{ 9 a0 (x) 

2UU.nl a-" 



+ nl-^Ar)&^B%g(x)+ (5.22) 



where [y k y l ] = y k y l — r 2 5 kl /3, and we assume that the third term with the tensor-type 
operator BjF^ appears first at t^f{> 0) loop of the perturbative expansion in the OPE. 
Therefore, with this assumption, we have 

$1 = _ e si + A T J ~ 24 (5 . 23) 
T T 2(33 -2iV» V ' 

where Aj? = 7r/(2d) with being the anomalous dimension of the operator B^ 9 ^. The 
calculation of anomalous dimensions for all 6-quark operators in the previous section shows 
that A T < 24, so that /3f 7 < ^ < 0. 

An extra matrix element we need is given as 



[y k y l ] /n 

XCT— — (0 



3 D 1 ,l) = B T( (>( 3 D 1 ) , (5.24) 



where Bj- is a further non-perturbative constant. 

Using the above results, we have for (S,I) = (1,0) 



fi =1 (y) * {B VI + (-logrfo°B II }<p( 3 S 1 )+(-logrfTB T 4>{ 3 D 1 ) . (5.25) 
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By applying V 2 we obtain 



2m raw r z v ' 



From V(y)y J £ l {y) = (E + V 2 /(2m))(^ =1 (j/), we obtain 



6H^^ ( s^). (5 . 26) 



y c 10 (r)^F 10 (r)^, (5.27) 
n vi 

Vr(r)^^°(r)^-, (5.28) 



where 



F-(r) = ^^p. (5.29) 

This shows that the central potential V c (r) diverges as F w (r) in the r — > limit, which is 
a little weaker than 1/r 2 , while the tensor potential Vr{r) diverges as F^°(r) in this limit, 
which is not stronger than F 10 (r) since /3q° — 1 > /3<r°- 

6. Evaluation of matrix elements 

We rewrite 3-quark operators in terms of left- and right- handed component: 

B f x , a = Bif;(C l5 P x )^(zT 2 ) gh , (6.1) 

B f XY,a = (Px)a B^, (6.2) 

for X, Y = R or L. In terms of these we have 

B{B 9 B = [Bj + Bjj + Bjjj + B IV + B V + B VI ] f J 8 , (6.3) 



where 



= [BrrBrr + B LL B LL \% , (6.4) 

( B n)% = [BrrBrl + B RL B RR + B LL B LR + B LR B LL \% , (6.5) 

i. B ni)% = [BrrB L r + B L rBrr + + B RL B LL \% , (6.6) 

= [Brl Brl + B LR B LR \% , (6.7) 

(B v )% = [Br L B L r + B LR B RL \% , (6.8) 

(^/)^ = + B ll B R r\% . (6.9) 



Note that we take (x, i) = (0, 0) in the above operators. We need to know 

> a/3 



(0\(Bi) f a s B \2N,E) (6.10) 
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for i = II, VI. 

For f ^ g, Lorentz covariance leads to 



(0\B f x B°r\2N,E) = Yl C x §(s)P A u(p,a 1 )P B u(-p,a 2 ) 



.11) 



A,B=R,L 



where s = E 2 = ^yp 2, + m N with the total energy E in the center of mass frame, Oi 
(i = 1,2) is the spin of the i-th nucleon, and C XY is an unknown function of s. Using 
invariance of QCD under the parity transformation PBxP~ l = 1aB x where R = L and 
L = R, we rewrite eq.( 6.11 ) as 

Q = (U\PB x B 9 Y P- l P\2-H,E) =J2Cx?PajM-P,*i)PbJMp,<72) 

A.B 

= E C M P MP, °i)Pbu{-p, a 2 ), (6.12) 



A,B 



where 74it(— p, <7i) = u{p,a{) is used. The above relation implies Cj^ = C XY . Using this 
property for the unknown functions C XY , we have 



<°l (Bn)lf 9f \2N, E) = C*£f LR {(P R ®P R + P L ® P L )u(p, a^-p, a 2 )} a p T(3a (6.13) 



and 



(0\(Bvi) f a f 9f \2N,E) = C^UPr ®P L + P L ® P R )u(p,a x )u{-p,a 2 )} a ^ a (6.14) 
Taking p = (0,0, p z > 0) and Dirac representation for 7 matrices Q, we have 



u(±p, +) 



1 



/ E N + m N \ 




^JE N + niN 



V 



TPz 





u(±p, 



1 



y/E N + m N 



J 



( \ 

En + mN 




(6.15) 



where En = y p 1 + m.^. For 7 = 1 ( fg + gf) and S = (eri = + and 02 
explicit form for the spinors gives 

{(P R ®P R + P L ® P L )u(p, +)u{-p, -)}i 2 -2i = E N , 
{(P R ®P L + P L ® P R )u(p, +)u(-p, -)}l2-21 = m N , 

while, for I = ( fg — 5/) and S 1 = 1 (07 = + and a 2 = + ), we have 

{(P fl ® P R + P L ® P L )u(p, +)u(-p, +)}n = miv , 
{(P fl ®P L + P L ® P R )u(p, +)u(-p, +)}n = Ejv • 



We finally obtain 



(0\(B VI ){ 9 2 +9f \2N,E) rn N C ^(s 



RL, 



the above 

(6.16) 
(6.17) 



3.18) 
3.19) 



(6.20) 
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for fg + gf and ((Ti,a 2 ) = (+,-) ( ), and 

(0\(B H ){(- fa \2N,E) _m N C^- LR {s) 



(6.21) 



(0\(B VI )(r f9 \2N,E) E N C^-(s) 

for fg - gf and (ai,a 2 ) = (+, +) ( 3 Si ), where s = 4E%. 

Unfortunately, we can not determine the sign of the ratio for these matrix elements. As 
a very crude estimation, we consider the non-relativistic expansion for constituent quarks 
whose mass toq is given by uiq = mjv/3. In the large tuq limit, 74^0 = Qo and 741*0 = no, 
where a subscript for q and u means the 0-th order in the non-relativistic expansion. 
In this limit, it is easy to show Cj& = C for all X,Y,A,B, so that C^ +LR = 2C and 
Cjj^ = C . Furthermore the first order correction to C^y = C vanishes in the expansion. 
Therefore in the leading order of the non-relativistic expansion, we have 



(0\(B VI ){l +9f \2N,E) 
for (0-1,02) = (+,-) ( 1 S'o ), and 

(0\(B n ){(- f9 \2N,E) 



(0|(fr/)(f +g/ |2N,S) a2 + Q (i) (6 22) 



(0\(B VI ){°- fg \2N,E) 




(6.23) 



for (01,(72) = (+, +) ( 3 5*i ). For both cases, we have positive sign for the ratio, which 
gives repulsion at short distance, the repulsive core. 

7. Conclusions and discussion 

The OPE analysis leads to conclusion that the S-state potential at short distance behaves 
as in (|l.3| ) with (|l.4| ). However perturbative considerations alone can not tell the crucial 
sign of the overall coefficient Ce- Moreover we found that the latter was also not directly 
predicted by chiral PT . A crude estimation using non-relativistic quarks suggests that Ce 
is positive, hence predicting a repulsive core, which diverges a little weaker than r~ 2 at 
small r. The leading corrections involve small powers of logs and hence it could happen 
that the dominant asymptotic behavior appears only at extremely short distances. 

Our analysis suggests that the repulsion of the NN potential at short distance is related 
to the difference of anomalous dimensions between a 6-quark operator and two 3-quark 
operators at 1-loop, and to the structure of the composite operators which probe the NN 
states. The explicit 1-loop calculation indicates that a combination of fermi statistics 
for quarks and the particular structure of the one gluon exchange interaction determines 
the sign and size of the /3's. The appearance of zero effective gamma eigenvalues is simply 
explained by chiral symmetry, however we were unable to find a simple proof of the absence 
of positive eigenvalues established by explicit calculation. One speculation is that the latter 
intriguing pattern can be explained by the relation between the 1-loop QCD anomalous 
dimensions and those of super YM theory [10]. 
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At higher order in the perturbative expansion, tensor operators appear in the OPE. 
Using this fact, we also found that the tensor potential also diverges a little weaker than 
r~ 2 as r — > 0. 

There are several interesting extensions of the analysis using the OPE. An application 
to the 3-flavor case may reveal the nature of the repulsive core in the baryon-baryon 
potentials. Since quark masses can be neglected in our OPE analysis, the calculation 
can be done in the exact SU(3) symmetric limit. It is also interesting to investigate the 
existence or the absence of the repulsive core in the 3-body nucleon potential. Such an 
investigation would require the calculation of anomalous dimensions of 9-quark operators 
at 2-loop level. Certainly more precise evaluations (also involving numerical simulations) 
of matrix elements {0\Ox\E) will also be needed to theoretically predict the nature of the 
core of the NN potential. 
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A. Explicit calculations of the divergent part for 6— quark operators at 
1-loop 

In determining the divergent parts for 6-quark operators we will consider the various cases 
in turn adopting the following mechanical (if rather inelegant) procedure. We first list the 
various operators which can appear and determine their linear relations. Then we compute 
the divergent parts using ( |4,25| ), initially keeping iV explicit in the formulae to indicate 
from which part of ( [4.25| ) the terms originate. Finally we set N = 3 and use the constraint 
equations to express the result in terms of linearly independent operators. 



A.l S = and 1=1 case 

1. Bf = B f {° ,B f J° . + B f J? R f Ji , 

I a[/3a] p\pa] a[/3a] $]/}&] 

We prepare following 11 operators 



B 5 = B'JZBft', B 6 = B^B%, B 7 = B^B%% B 8 = B^B^, 

B9 = B^BlZ, B W = B f J/XZ> B ^ = B L% B lZ . (A- 1 ) 

which are all the possible operators corresponding to the Dirac labels 111222 in Table || 
(in the 4/2g case). In terms of these we write 

Bf J[Sa) B i{L\ =B 1 +B 2 -B 3 -B 4 . (A.2) 
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Taking -Bi,2,3,4 as independent operators, the constraints from ( 4.44Q are: 

-3Si 



B 5 
B 9 



B4 — 4B3, Bq = B\, B 7 = 2B 3 — B4, B$ 
3B 2 , Bio = B4 — 2B 3 , Bn = —B 2 ■ 



(A.3) 



The 1-loop corrections to -61,2,3,4 can be calculated from To as 
rj 1 ) = (5 - 9N)Bi - 4NB 3 - 2{N + 1)B 4 - 2B 5 + ^(-3^6 - 2B 7 + B 8 + 4B W ) 



2 
3 



(5 - 9N)B 2 - 4NB 3 ~ 2 ( N + 1 ) 5 4 - 2jB 5 + N ( 2B io + 3B n - 4B 7 - B 9 ) , 
(1 - UN)B a - 2N(Bi + B 2 ) + (2 - N){B 4 - B 5 ) + 2N(B W + B n - B 6 - B 7 ) , 
(3 - 8N)B± - 4(JV + l){Bi + B 2 ) + 4(2 - N)B 3 + (N - 2)B 5 + 2N(B W - B 7 ) , 

(A.4) 



where an overall factor d/e is dropped on the rhs for simplicity, i.e. using the gauge 
identities eq. ( |A.3j ) we obtain 

(A.5) 



r, (1) = jlii B f , 



with the matrix 7 given by 



So we obtain 



7 /(2d) 



/-40 -40 8 \ 
-40 -40 8 
-12 -12 -60 12 
\-16 -16 -32 -8 J 



(A.6) 



1+2-3-4 



12~-Bi + 2-3-4 



(A.7) 



For B 



ffg gffg 



the result can be obtained from the above by the interchange of a and 



a. In terms of Bi = Bi(a •<-)• a) , we have 



rf) w = (b + b) 



1+2-3-4 



1— loop,div 



-12-(£ + 5)1+2-3-4 



(A.8) 



This corresponds to the result for the first 1 = 1 operator for the entry 111222 in Table |2[ 
The other eigenvalues of 7 are obtained similarly. 



2. BO} 



B S \ 9 R ,B S \? , + B n J ,B J J® , + (a,P++&, $) 



iffg r>ff 9 



We prepare the following 6 operators, 



Bi = B f J R 9 B f Jl „ B 2 = B ffa f,B f Ji „ B 3 = B J R J "B 



Ba 



Tiff 9 rtffg r> _ r>fff r>9f9 



in terms of which we have 



ffg r>ffg 

0P<* a[$&] 



r _ Tjfff r>gfg 



(A.9) 



B ffg ,,ffg . i>fJg n ff-> 

a\Ba\ ft flA 



a[/3&] PIM 



Bi — B 2 — B 3 + i?4 . 



(A.10) 
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The 1-loop corrections to #1,2,3,4 can be calculated as 



1 



(3 - N)B 1 - 2(N + l)B 2 + (N — 2)B 3 - 2B 4 + NB 5 - 2NB 6 , 
AB 2 - 4(iV + 1)B 1 + 2(iV - 2)B 4 - NB 6 , 



rf } = AB 3 - 4(N + l)B 4 + 2(N - 2)Si - NB 5 , 

rf } = (3 - N)B 4 - 2(N + 1)5 3 + (AT - 2)B 2 - 2B 1 + NB 6 - 2NB 5 . 

Therefore we have 

r(l) -19 d R 

1 1-2-3+4 — i2 -^l-2-3+4 • 

For B,~ M*J 9 , + ,B{{~ ,, we introduce -Bj = -Bj(a ++ a), so that we have 

a [/3a] /3[/3a] a[fia] p[p&]> V ; ' 



= (-B + -B) i_2-3+4 



1— loop,div 



12-(B + B)i_2_ 3 +4. 



(A.11) 
(A.12) 

(A.13) 



We prepare the following 15 operators: 



Si = B ff9 B f J 9 , B 2 

1 rvnrv qq~. * 



B ff9 B f J 9 , B 3 



B ff9 B ff9 B = B ff9 jiff 9 

aa P /3 Pa' 4 



R Bl jy , B 5 

aa P paP' 



Tjffg nffg 
n pPa n p aa ' 



B « = Bf 4l B & B ? = b ^ b Z 9 p' c - = c * 

Dl = B laP B iZ, D * 



B fgg B fff ,) = B fgg B fff D 

Paa PolP 6 a PP Paa> 4 



_ r>fgg r>fff 

- n aap n p a/3 

in terms of which we write 

B f { 9 M f J 9 , = Bi - B 2 - B 3 + B 4 . 

p[p a ] 1 2 

Taking #1,2,3,4,5 as independent operators, others can be expressed as 



jjfgg rjfff j~) _ rjfgg Tjfff 



(A.14) 



(A.15) 



Bq = —B 2 , B 7 = B 3 — 4E?2, 

C\ = —2B\ — i?5, C 2 = D\ = —D 4 = 2B 2 — B 3 , C 3 = —3B 4 , 
D 2 = —B^,D 3 = —B 4 , D5 = B5 — 2B\. 

The 1-loop corrections to -61,2,3,4 can be calculated as 

rj 1} = -±NB 1 - 2(N + 1){B 2 + B 3 ) + (N- 2)B 5 - 2(B 6 + B 7 ) 
+ N(d - 2C 2 - 2D 1 -D 2 + D 5 ) 
= -24Bi - 22B 2 + 25 3 + 4B 5 , 

= -5NB 2 - 2(N + l)(£i + B 4 ) - 2B 5 + (N - 2){B % + B 7 ) 

+ N(D 4 -C 2 -D l - 2D 2 ) 

= -8B 1 - 38B 2 + 10B 3 - 8B 4 + 4B 5 , 

(2 - 4A0£ 3 - 2(7V + 1)(B 4 + 2B ± ) + (N — 2)(2B 6 + B 7 ) - N(C 2 + 2D L ) 
-I6B1 - 24£ 2 - 8B 4 , 

(2 - N)B 4 - 2{N + 1)(B 3 + 2B 2 ) - 2{2B 6 + B 7 ) + N(C 3 - 2C 2 + 2D 3 - 4£>i) 

(A.16) 



-40B 2 + 8B 3 - 16B 4 
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We therefore obtain 



r (i) 

1 1-2-3+4 



We have 3 more structures: B 



ffg B ffg 



D 



0. 

ffg B ffg 



and B 



ffg B ffg 



(A.17) 
We introduce 



a[fSa] /3[/3a]' &[$&] p[$a] a[0a] $[$& 

operators B[ = —Bi(a o /3) for the first one, Bi = Bi(u a) for the second, and 
= i?j'(a -H- a) for the third. We then have 



(r?y (1) = o. 



(A.18) 



4 R 01 



_g//s _|_ B ffg B fjg 



a{(3a] 0[/3a] 

We prepare the following 12 operators: 



/,', = i! ff " ll f f:l . B 2 = B ff £B f /A, B 3 



a/3a 13a 1 



a/3a /3A/3' 

— B^J~B^ 99 



b a/3 /3 Z 3 ""' 1 a/3/3 Z 3 ""' 2 

D 2 = B f9 AB f U, D 3 = B f ".B ff l, 



B ffg B ffg B 

aa/3 1313a 1 4 

B ff A B f99 . (\ 

aaa aaa' W J 



: B ffg . B ff9 , B 5 = B ff AB ff A, 

aa/3 /3a/3' b aaa /3/3/3' 

B fff B fgg D B fgg B fff 

aa/3 Pa/3' a/3/3 P aa ' 



in terms of which we write 



B ffg B ffg 



— B\ — B 2 — B 3 + i?4 . 



(A.19) 



a[/3a] P\p&] 

The constraint from gauge invariance leads to 

AB X + d + D 2 = 0, 4B 4 + C 2 + D 1 = 0, B 2 + B 3 + B 5 + B 6 + C 3 + D 3 = 0.(A.20) 
The 1-loop corrections can be calculated as 



p(l) 

1 1 


= (5 


- 3N)B 1 


- 2(N + 1){B 2 + B 3 ) - 2(£? 5 + -B 6 ) + N(D, 


, _ 2D 3 + Ci - 2C 3 ) , 


p(l) 

1 2 


= (3 


- AN)B 2 


- 2(N + l)(Bi + B 4 ) + (JV - 2)(B 5 + £ 6 - 


B 3 ) - N(C 3 + D3) , 


r (l) 
1 3 


= (3 


- AN)B 3 


- 2(iV + + B 4 ) + (iV - 2)(S 5 + B 6 - 


B 2 ) - N(C 3 + L> 3 ) , 


r (l) 
1 4 


= (5 


- W)B 4 


- 2(iV + 1)(B 2 + B 3 ) - 2(B 5 + B 6 ) + iV(D] 


-2D 3 + C 2 - 2C 3 ) . 



Therefore we have 



.(i) 

1-2-3+4 

ffg rjffg 



= 0. 



We introduce B~ = BAa. <+>■ a) for Bi!x ,B'J 9 so that we have 



01 



0. 



(A.21) 



(A.22) 



(A.23) 



5. flOi 



o//fl r//s r//S o//s „„j r01 
a[/3&] /3[/?a]^<*[/H^?[/Sfi] ^VJ 



r//9 r//3 r//S r//9 
^a^a] D p[pa\ a[$a] W<A 
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In this case we have to prepare 29 operators as follows: 



Bi 


afia /§/9a' 


B 2 


— B^ 9 B^ 9 

apa PaP" 1 


B 3 


— B^ 9 B^ 9 

aa/3 PPa 1 


B 4 


— gffg gffg 

a&P Pap 1 


B 5 


— B^ 9 B^ 9 

PPa paa'' 




0a$ Paa'' 


B 7 


= B ff AB f J 9 , 

aPP a P a ' 


B 8 


= B f l 9 B f J 9 , 

aPP aa P' 


Ci 


= B ff AB f99 , 

aPP P«"' 


c 2 


= B ff lB f99 

aPP aa P' 


c 3 


— B^ 99 

aa/3 /3/9a' 


c 4 


_ gfffgfgg 

aap a$P' 


Di 


— ftfgg gfff 

Paa Pap 1 


D 2 


_ ftfgg gfff 

PPa Paa' 


D 3 


_ gfgg gfff 

aap paP 




— B^ 99 B^^ 

&P$ Pact' 




— r//9 r//9 


x 2 


— B^ 9 B^ 9 

a P a PPa 


, x 3 


__ ufgg r>fff 

aaa pppi 


Yi 


_ off 9 -off 9 
~ °ppa 


Y 2 


— Rff9 £}ff9 

aaa pfap 


Zi 


= B ff AB ff J, 

&PP a/ 3 "' 


z 2 


— gffg gffg 

aPP aPa 1 


Z3-- 


- pjffg r>ffg 

aaa ppp' 


Z 4 ~- 


= B ff J R Bl 99 , 

aa P PPa' 


z h 


— B^ 99 B^^ 

aap PPa 1 


Vi 


— R^'R^ 9 

&PP "a/ 3 ' 


v 2 -- 


- B^ 9 B^ 9 

aaf3 &PP' 


V 3 = 


_ f>fffr>fgg 

aaa pfcp 







in terms of which we have 

=Br~B 2 -B, + B^ B 5 , (A.24) 

B ^PM B m A] = +Y 1 + Z 1 -V 1 = B 9 . (A.25) 

The constraint from gauge invariance gives 

C\ = —2B\ — Xi, D\ = -B\ — B 5 — Xi, X 2 = —2(B\ + B 5 ) - X\, X 3 = 2B 5 + Xi, 
D 2 = C 2 = -2B 2 - Y lt Y 2 = -4B 2 - Y lt Z 3 = Z 2 + B 6 - C 3 , Z 4 = D 3 + Z l - Z 2 , 
Z 5 = B 6 + Z 1 - Z 3 , B 6 + D 3 = B 7 + C 3 , B 3 + Z 1 + B 6 + D 3 = 0, 
C 4 = -B 4 -B 8 - Vy, D 4 = -2B 4 - V u V 2 = -2B 4 - 2B 8 - Vy, V 3 = 2B 8 + Vi,(A.26) 

where By : 2,3,4,5,6, 7,8) Xy, Yy, Zy j2 and V\ are taken to be independent. The 1-loop corrections 
to -£^2,3,4 can be calculated as 

rf } = (2 - N)B X - 2(N + l)(B 2 + B 3 ) + (N — 2)B 5 - 2B 7 + N(d - 2C 3 - D 1 ) , 

= 2B 2 - 2(N + + B 4 ) - 2{B b + B 8 ) + N(C 2 - 2C 4 + D 2 - 2D ± ) , 

rf } = (2 - 2N)B 3 - 2(N + + B 4 ) + (N - 2)(B 6 + B 7 ) - N(C 3 + D 3 ) , 

V{ 1) = (2 - N)B 4 - 2{N + 1)(B 2 + B 3 ) - 2B 6 + (N - 2)B 8 + N(-C 4 - 2D 3 + L> 4 ) , 

(A.27) 

which leads to 

r i-2-3+4 = 3(iV + 2)(5i + B 4 ) - (AN + 6)B 2 - (2N + 6)B 3 + N(B 5 - B 6 - B 7 + B 8 ) 
+ AT(Ci -C 2 -C 3 + C 4 + D 4 - D 2 - D 3 + D 4 ) . (A.28) 

On the other hand, constraints give 

d - C 2 - C 3 + C 4 + Dx - D 2 - D 3 + D 4 = -3Bi + AB 2 + 2B 3 - 3B 4 

- B 5 + B 6 + B 7 - B 8 + 2(-X 1 +Y 1 + Z 1 - V x ) , (A.29) 

so we finally obtain 

ri 1) 2 -3+4 = 6^65 + As % . (A.30) 
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We introduce operators Bi = Bi(a -H- a) for B^ : q, so that we have 



(r^) (1) =6-(S + Z?) 5 + 6-(£ + Z?) 6 , (r° y 1 / ) (1) = 24-(S + fi) 6 . (A.31) 



A.2 5 



1 and 7 = case 




a\p,a] a[8,a] &\J3,&] &\J3,&] 



We prepare the following 11 operators: 




C2 = BfJ a iB%% C 3 



r>f.fg ptggf 
p>fff p>ggg 



in terms of which 




B\ — B2 — B 3 + i?4 . 



(A.32) 



Constraints are given by 



C\ — C2 — 2B\ — B4, C5 — B4 — 4Bi, C3 — B2, Cq — SB2, C4 — B 3 , C-j — — 3B3. 
tie 1-loop corrections can be calculated as 

if } = -(3 + 7N)B l - 2(B 2 + B 3 ) + N(B 4 - C 5 ) - 2(d + C 2 + C 3 + C 4 ) , 
if } = -(1 + 3^)52 - 4Bi - 2(N + 1)B 4 - 2Ci - AC 2 - (N + 2)C 3 + iV(C 6 - 2C 5 ) , 
if } = -(1 + 3iV)B 3 - 4Bi - 2(AT + i)5 4 _ 4C1 - 2C 2 - (AT + 2)C 4 + A^(C 7 - 2C 5 ) , 
rf } = (1 - 6N)B 4 - 4(AT + i)(B 2 + £ 3 ) + 4NBi + 2(AT _ 2)(Ci + C 2 ) - NC 5 , (A.33) 



For B,? ,B",{ the result can be obtained from the above by the interchange of a and 

a[$,a] &[$,&]' 

a. In terms of B^ = B^a ■<->■ a), we have 



and therefore we have 




— — 4-f?i_2-3+4 . 



e 



(A.34) 



// a[/3a] of/Ja] d[/3a] a [/3a] V >K 



(A.35) 



We prepare the following 6 operators: 




-29- 



in terms of which 



B f f 9 M" r l , + B f {? B 99 { , = B 1 - B 2 + B 3 - B A . (A.36) 

a[M a [pa] a[pa] a[/3a] 1 2 4 V ; 



The 1-loop corrections can be calculated as 

if } = (AT + l)(Bi - 2B 2 ) - N(2B 3 + B 4 ) + iVC 2 - 2d , 

if } = (2JV + 3)B 2 - 4(A + l)B 1 - 2NB 3 + (N - 2)d , 

if } = (N + 1)(B 3 - 25 4 ) - N{2B± + B 2 ) + ACi - 2C 2 , 

if } = (27V + 3)B 4 - 4(A + 1)S 3 - 2iV5i + (iV - 2)C 2 , (A.37) 

and therefore we have 

r i-2+3-4 = 20^-2+3-4 • (A.38) 

For B ll$a] B9 4h] + B *[M B % & y we introduce B< =^(aHa), so that we have 

(r]?) {1) = 20^(5 + 5) 1 _ 2+3 _4. (A.39) 
3- = + + (a, /3 o a, /3) 

We prepare the following 15 operators: 

C 2 = B ff9 B? 9f , C = B f / 9 B f J 9 , C- = B"l B"l, C 5 = B" f R B ff9 , C* 6 = B"/ R B ff9 , 

z app acta' J /3pa tma' 4 apa apa^ aap apa' D app acta' 

C 7 = B^B 9 & %, C 8 = B^BIZ, C 9 = B 99 i a B%% C 10 = B^B 9 ? 9 ^ C u = BffljB™, 
in terms of which we write 

^& B % a] =B 1 -B 2 -B 3 + B 4 . (A.40) 

Constraints are given by 

2C 2 = 2B\ — B4 + Ci, 2C4 = —C\ — B 4l Cg = -B4 — 2i?i, C 7 = —C\ — 2B±, 
C5 = — C3, 6*8 = —2B 2 — C3, Cg = —2B 2 + C3, C10 = — 3i?3, Cn = — -B3. 

The 1-loop corrections can be calculated as 

rf } = -5NBi - 2(N + 1)(B 2 + S 3 ) + (N - 2)(Ci + C 2 ) - 2C 3 

+ A(-C 4 - 2C 5 + C 6 - C 7 ) , 
rf } = -4AB 2 - 2(A + l)(Bi + B 4 ) - 2(d + C 2 ) + (A - 2)C 3 

+ A(-2C 4 - C 5 - 2C 7 + C 8 + C 9 ) , 
if } = (2 - N)B 3 - 2(N + 1)(B 4 + 2Bi) - 2(Ci + 2C 2 ) + iV(-4C 4 - 2C 7 + Ci + 2C U ) , 
r { l ] = (2 - 4A0£ 4 " 2(A + 1)(B 3 + 2S 2 ) + (N - 2)(Ci + 2C 2 ) - iV(2C 4 + C 7 ) . (A.41) 
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We therefore obtain 



r (i) 

1 1-2-3+4 



0. 



(A.42) 



We have 3 more structures: BU,S ^B gg ,{ B^? B gg l , and B^ g B gg ( , . We introduce 

a[/3a] a[/3a}> a[/3a] a[/3&] a[/3a] 

operators B[ = —Bi(a «-»■ (3,f -H- g) for the first class, B{ = B{(a <H> a) for the second, and 
Bi = B'(a f-> a) for the 3rd. We then have 



10 \(i) 



(r&) 



o. 



4 R 10 



a[f3a] a[/3a] ' a[(3a] a[/3a] 

We prepare the following 16 operators: 



(A.43) 



n. = B ffg B gg l, ii, 

apa apa 



Co 



nffg r>ggf 

a/3/3 «aa' 3 



B ffg B gg f B 3 

apa aap 

B"I B ffg C 4 

a/3a apa 



Q 7 _ R SS/ r//9 

D 4 = B"lB ffg , 

4 a/9/3 aaa' 

in terms of which we write 



1313a aaal 8 



#/// £>ggg 



a/3 (3 1 



B ffg B ggf Bi 

aa/3 af3a 

gfffgggg q . 

aaf3 aa/3' 



D, = 



E}ffg £>ggf 

a/3a a(3a ' 



Do 



joffg r>ggf q _ r>ffg r>ggf 

aa(3 aa$' 1 a&& a/3/3' 



T^ggf rjffg r< 

B aaa B p^ U 



j^fffj^ggg 

af3(3 



a/3/3' 



af3a af3a 



fiff? £>ggf 

a[/3a] a[/3a] 



= B\ — B2 — B3 + B4 . 



Constraints give 

d + C 2 + C 3 + C 4 = -(Si + £ 4 ), C 5 + C 6 + C 7 + C 8 = -4(B 2 + £3) 
The 1-loop corrections can be calculated as 



(A.44) 



(A.45) 



1 

2 



(1 - 2N)Bt - 2(N + 1)(B 2 + £3) + ^(£4 - C 3 - C 4 ) + (iV - 2)(d + C 2 ) , 
(1 + iV)(£ 2 - 2Bi - 25 4 ) + W(C 5 + C 6 - 2C 3 - 2C 4 ) - 2(Ci + C 2 ) , 
(1 + iY)(£ 3 - 2Bi - 2B 4 ) + iV(C 7 + C 8 - 2C 3 - 2C 4 ) - 2(d + C 2 ) , 



= (1 - 2iV)5 4 - 2(7V + 1)(B 2 + £3) + N(Bi - C 3 - C 4 ) + (iV - 2)(d + C 2 ) , 

(A.46) 

which give 



1-2-3+4 



8-B 



1-2-3+4 • 



(A.47) 



We introduce B, t = BAa <+>■ a) for £(f„ 9 so that 

' 'V ' a[/3a] a[/3a]' 



10 \ (!) 



8-(5 + 5) 1-2-3+4 • 



(A.48) 
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5. B$ 



B 



Jfg oggf , off 9 oggf QTir i Rio _ off 9 oggf , offg oggf 



iff 9 oggf 



We prepare the following 30 operators: 



Bi 


= B ff /B 99 {, B 2 

a p & a/3a> * 


= B ff9 B? 9f R , 

a/3a ««P' 


B 3 


= B ff9 ~B 99 {, B 4 = 

aa/3 a l 3a ' 4 


B ffg B ggf 

aa/3 


B 5 = 


B ff9 B 99f 

«P a d/3d' 


B 6 


= B f J 9 B" f ~, B 7 

a P a aa/3 1 ' 


= B ff9 B"J , 

aa P a/3a' 


B 8 


= B f J 9 B"{, Ci = 


Qggf gffg 

0&0 


c 2 = 


B" f .B ff9 , 


c 3 


_ oggf offg r< 


= B™fB ff9 

aj3a 


c 5 


= B f ( 9 B 99f , Cq = 


offg oggf 

««« /3/3a' 


c 7 = 


offg oggf 

0a0 


c 8 


— offg oggf c 1 


_ oggf offg 


C10 


._ oggf offg /"» . 


_ oggf offg 

aaa 00a 


j C12 






= B fff .B 99 A, D 2 

aa/3 P aa ' z 


= B 99 J B f l 9 , 

f3aa a p & > 


D 3 


= B ff9 ~B"l, I), = 

a/3/3 aaal 


. gffg £>ggf 

aaa 


D 5 = 


--B"/B f R f A, 

af3a POQ' 


#6 


= B f J/-,B 999 D 7 


= B ff9 B? 9 J R , 

a/3a ««P' 


D 8 


= B ff (B 999 ~ , Da = 

a pp aaa> » 


-. B ff9 ~B 99 l, 

aa/3 P« Q ' 


D w 


_ gfff gggg 

aaa p & 0' 



in terms of which we have 



a[/3d] + a a[0a] a &[ p &] - » 5 + C 6 , 

£>5 = -Bl-2-3+4 , £>6 = -B5-6-7+8 • 



(A.49) 

(A.50) 



Constraints give 



C\ = C3 = B% — 2Bi, D\ = B$ — 4Bi, D 2 = —B\, C 7 + C% + Cg + C\o = 2B^ — 2B4, 
C 2 + C 5 + C11 = 2B 6 - 3B 2 , C 4 + C 6 + C12 = 2,67 - 3^3. (A.51) 

The 1-loop corrections can be calculated as 



- 2B l - 2{N + 1)(B 2 + B 3 ) - 2{C 5 + C 6 ) + N(d + C 3 - 2C 2 - 2C A ) , 

2 

r (i) 

1 3 



T^' = (2 - iV)B 2 - 2(iV + l)(5i + B 4 ) + {N- 2)C 5 - 2C 8 + N{-C 2 + Cn - 2C W ) 



(2 - JV)B 3 - 2(N + + B 4 ) + (N- 2)C 6 - 2C 7 + N(-C 4 + C12 - 2C 9 ) , 

if } = (2 - 2^)54 - 2(N + 1)(B 2 + S 3 ) + (N — 2){C 7 + C 8 ) - N(C 9 + C 10 ) , (A.52) 



which give 



(A.53) 



We introduce operators Bi = Bi(a <H> a) for ^5^, so that we have 

= 6-(B + B) 5 + 6-(B + B) 6 , (r^) (1) =24-(fi + £) 6 . (A.54) 
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Dirac indices 


lj /{2d) 


I 


mill 


-24 





111112 


-24 


0,1 


111122 


-4 







-24 


0,1 




-40 


2 


111222 


-4 







-12 


1 




-24 


0,1 




-40 


2 




-72 


3 


111113 


-16 


0,1 


111123 


-6 


0,1 




-16 


0,1 




-24 


1,2 



Table 2: Eigenvalues 7j of the anomalous dimension matrix 7 and isospins of the corresponding eigenvec- 
tors for the case 3f3g. 
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Dirac indices 


7i/(2d) 

< J 1 \ t 


/ 


111223 





0,1 




-6 


0,1 




-16 


0,1 




-18 


1,2 




-24 


1,2 




-48 


2,3 


111133 


-4 







-16 


0,1,2 


111134 





1 




-4 







-12 


1 




-16 


0,1,2 


111233 


4 


1 




-4 







-8 


0,1,1 




-16 


0,1,2 




-32 


1,2,3 


111234 


20 







8 


1 




4 


1 







1 




-4 







-8 


0,1,1,2 




-12 


1 




-16 


0,0,0,1,1,2,2,2 




-32 


1,2,3 


112233 


8 







4 


1 




-4 


0,0,1,2 




-8 


0,1,1,2 




-16 


0,1,2 




-28 


2 




-30 


1,2,3 



Table 3: 



As in Table 
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2 (continued). 



uirac maices 




J 
1 


llzzo4 


on 
zU 


n 
U 




1 o 
Iz 

Q 

o 


1 

U, 1 




4 


i 




n 
U 


i> 1 




-4 


0,0, 1,2 




Q 

— O 


n i 1 9 

0, 1, 1, 1 




-12 


1 




i fi 

— ID 


n n i 1 999 

U, U, 1, 1, z, z, z 




— Zo 


o 
z 




39 
— oz 


19 3 

1, Z, ,5 
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1,2,3 


lllooo 


— D 


n i 

U, 1 




94 
— Z^ 


n 1 9 ^ 

U, 1, Z, 6 


111 'X'XA 
111004 


n 
u 


n i i 9 

U, 1, 1, z 




— 


n i 

0, 1 




1 8 


1 9 
1, 2 




OA 
— Z4 


n 1 9 ^ 

U, 1, Z, cS 


1 1 9QQ/I 


9/( 
Z4 


n i 

0, 1 




6 


0,1 







0,0,1,1,1,1,2,2 




-6 


0,1 




-12 


1,1,2,2 




-18 


1,1,2,2 




-24 


0,1,2,3 




-30 


0,1,2,3 



Table 4: As in Table 2 (continued). 
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B. Some useful formulae for angular momentum states 



B.l Eigenstates 

At given J, there are 2 distinct states, the spin-singlet (5 = 0) state and the spin-triplet 
(5 = 1) state. 

The singlet state is denoted as 1 Jj, since it has 5 = and J = L. The fact that 
I + L +- S must be odd to satisfy fermion anti-symmetry gives 1 = for odd J and 1 = 1 
for even J. The eigenstate with J z can be easily obtained as 



fjj.J*) = \ J g ,0)j, 



5=0 



(B.l) 



where \ J Z , S z )j,s = \ Jz)j <8> \S z )s- 

The spin-triplet state is classified into 3 types: 3 Jj, 3 (J± l)j. For the first one, 7 = 
(even J) or I = 1 (odd J), and vice versa for the other two types. By the Wigner-Eckart 
theorem, the matrix elements of the five operators do not depend on J z . Therefore it is 
enough to know eigenstates with J z = J only. Explicitly we have 



! J j, J) 



hi 



J -IX 



•Al 



VJ\J,0)j,i} 



\ 3 (J-1)j,J) = |J-1, !)./_!, 

1 



\J + I)j,j) 



:{|J-1,1) 



J +1,1 



^(J + l)(2J + 3) 

+ V2TT1 [v^7Ti)|J + 1, -1) - 1 J, o) J+ i,i] } . 

B.2 Evaluation of each operator 

Using these eigenstates, it is easy to see 



(B.2) 
(B.3) 

(B.4) 



<?i • <7 2 = 25(5 + 1) - 3 = -3, 1, 1, 1 

J(J+1)-L(L + 1)- 5(5+1) 



L-S 



0, -1, J-l, -(J + 2) 



(B.5) 
(B.6) 



for 1 Jj, 3 Jj, 3 (J — I) j and 3 (J + l)j, respectively. 

For 5i2 defined in (|5.2[) the results are more complicated due to the mixing between 
3 ( J — l)j and 3 (J + l)j. After a little algebra we obtain, 



/ 2(J-1) 6y/J(J+l) \ 



S 



12 



0, 2, 



2J + 1 



2J + 1 



6^/J(J + l) 2(J + 2) 



V 2J + 1 



(B.7) 



2J + 1 / 



C. The 1 = 2 2-pion system 



Here we consider the operator product expansion of two iso-vector pseudoscalar densities 
in QCD. 
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C.l Anomalous dimensions 

The local composite operators with 7r + 7r + quantum numbers in QCD with lowest dimension 

are 4-quark operators with dimension 6. There are 5 independent such (bare) scalar 
operators 

01 = dj^u ■ dj^u + d^jsu ■ drjfP/^u, (C.l) 

02 = du ■ du — d^u ■ d^^u, (C2) 

3 = Jy u • d-y^u - d7 M 75« • d'y^u, (C.3) 

04 = du ■ du + c?75U • d^^u, (C4) 

5 = do-^u • Jo>u, (C.5) 

where in this appendix we use the notation u = q u , d = q d and suppress explicit color and 
Dirac indices of the quark fields. 

There are also 3 independent such (bare) traceless tensors operators 

T» v = d^u ■ dj u u + dj>*j 5 u ■ dj u j 5 u - g^Oi , (C.6) 

= d^u ■ dj u u - dj^ 5 u ■ dYlsu - g^0 3 , (C.7) 
1 

D 

As in the main text, operators are renormalized according to the formula 

,2 



rf = da^u ■ da\u - - g^Os. (C.8) 



°A = ° A ~ 32^~e lAB ° B + ' ' ' ' (C,9) 



(R) _ nA 9 

and similarly for tensor fields. The results for the one- loop anomalous dimensions of the 
scalar fields can be found in [11]. The non- vanishing entries of the mixing matrix for the 
scalar case are 

7ii = -4, (CIO) 

722 = 16, (C.ll) 

733 = -2, 732 = -12, (C.12) 

744 = 10, 745 = 1/3, (C.13) 

754 = -20, 755 = -34/3. (C.14) 

We have extended the analysis of ref. p 1[ [ to the tensor case. Here we find 

7? = -8/3, (C.15) 

7^ = 2/3, 7? = "2, (C.16) 

7!? = -16/3. (C.17) 

It is useful to introduce the (onedoop) diagonally renormalized operators. For the scalar 
case we have 

2 

Xf ] =X A - -^iaXa + • • • , (C.18) 
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where 



Here 



*i = 0i, 71 = -4, (C.19) 

X 2 = 2 , 72 = 16, (C.20) 

X 3 = 3 + 2 -0 2 , 73 = -2 , (C.21) 

X A = e> 4 + B0 5 , 74 = (V964-2)/3, (C.22) 

X 5 = 5 + C0 4 , 75 = -(V964 + 2)/3. (C.23) 



B = — — y_^L = 0.01586, C = 32 - \/964 = 0.95165 . (C.24) 



30 

C.2 OPE for I = 2 vr-vr scattering 

The OPE for two ir fields can be written in QCD as 

n(x)7T(0) = J2la(x)B a + --- (C.25) 

a 

Here tt(x) is the field annihilating ir + , B a are the renormalized doubly charged dimen- 
sion 6 operators discussed in the previous section, x is spacelike (and for simplicity we 
assume its time component vanishes) and 7 Q (x) are c-number coefficient functions. [Here 
we "pretend" all operators are scalar, although in fact three of them are symmetric trace- 
less tensors. Taking into account their tensor structure however does not change any of 
our conclusions here.] The dots stand for higher dimensional operators with less singular 
coefficient functions. 

The short distance asymptotics of the 1 = 2 wave function is given by 

= (0|ir(x)7r(0)|2) + ( C - 26 ) 

a 

where 

B a = (0\B a \2) (C.27) 

are the (energy dependent) matrix elements of the local operators. 
In QCD we can write the divergence of the axial current as 

d^Afj, = (d^ l5 u) = m $o = m R <S> R , (C.28) 

where 

§ (x) = d(x) l5 u(x) (C.29) 

is a (bare) quark bilinear field with tt + quantum numbers and the quark mass parameter uiq 
is the sum of the u and d quark masses. [Here we used the fact that the axial current, being 
partially conserved, has renormalization constant Z = 1. There is a subtlety in dimensional 
regularization where because of the presence of 75 the renormalization constant is not equal 
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to unity. It is finite nevertheless and this does not alter our conclusions at the 1-loop level.] 
The (canonically normalized) pion field is defined as 

< x ) = -YT d »Mx) = ^r(x) , (C.30) 



n = ^. (c.3i) 



where Q is a constant: 

mlfn 

From this we see that the field <3?_r renormalizes with the inverse of the mass renormalization 
constant. 

The RG analysis of the pion-pion wave function goes along the same lines as in the main 
text for the nucleon-nucleon case. By inspecting the spectrum of anomalous dimensions 
we see that, again, only operators already present in the tree level expansion 

(Hu? = ~\o 2 + \o, = ~X 2 + ^-L— pQ - BX 5 ] (C.32) 

contribute to the leading short distance part of the wave function and even from this set 
we need only the operators with the largest anomalous dimensions. Since the coefficient of 
such an operator is asymptotically proportional to 

(-lnr)^ {A/A+do \ (C.33) 

where bo = 11 — 2Nf/3 and do = —16 comes from the mass renormalization. 
Numerically the spectrum of %4S is 

16 8 2 

(-4; 16; -2; 9.68; -11.02; -_;--;-), (C.34) 

corresponding to the spectrum of powers 

(-1.11; 0; -1; -0.35; -1.50; -1.19; -1.04; -0.85) (C.35) 

numerically. (Here we took iVf = 3 for simplicity.) 

We have, again, a leading zero eigenvalue and all the other powers are subleading. The 
next one is —0.35 so the wave function is asymptotically 

~ + ^ii~ b + ■■■ , (C.36) 

where b = 0.35. This corresponds to 
Here 

^ o = -^(O|0f)|2) (C.38) 
and ipi is proportional to (with a positive coefficient) the linear combination 

(0\0{ R) \2) +B(0\O { 5 R) \2). (C.39) 
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Note that the ratio ipi/il>o may be energy dependent. We need to calculate this ratio (or 
at least its sign) nonperturbatively, to be able to determine whether the potential in this 
channel is attractive or repulsive. ChPT is not applicable for this problem since there 
are too many extra low energy constants characterizing the matrix elements of 4-quark 
operators and in the end the sign of this ratio is left undetermined. In the absence of a 
reliable non-perturbative method to calculate the above matrix elements we try to estimate 
them by inserting a complete set of states in the middle of the operator and truncating the 
sum after the 1-particle contribution. This is very similar in spirit to the vacuum insertion 
method p2| , (oft rightly criticized) however surprisingly successfully applied to AS = 2 
weak matrix elements in the past. In this approximation 

(Q\(dT lU ) ■ (dT 2 u)\2) « (0|(dTm)|l)(l|(Jr 2 u)|2) (C.40) 

and therefore we have (in this approximation) 

(0\0[ R) \2) ^ -(0\O ( 2 R) \2) (C.41) 

and 

(0\O { 5 R) \2) ^0. (C.42) 

Thus the ratio ipi/ipo is positive in this naive approximation and the potential is repulsive, 
as indicated by the (quenched) lattice measurements |l3| . 
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